Two-electron correlated motion due to Coulomb repulsion 
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A Hubbard-type model is derived from the microscopic Schrodinger equation. We found that 
additional terms describing direct two-electron transitions must be added to the standard Hubbard 
Hamiltonian. Such a Hamiltonian generates two-electron pairing due to on-site Coulomb repulsion. 
We demonstrate that the electron pairs with opposite spin propagate across periodic structures 
via direct and sequential two-electron tunneling. This mechanism can be used for a generation of 
entangled electron current. Numerical calculations show stability of the electron pairs. 



PACS: 73.63.K, 03.65.X, 71.10.F 

Within the framework of the Hubbard model it has 
been earlier shown that the on-side electron-electron in- 
teraction (U) might generate coherent propagation of the 
two electrons 0,^. However, in the Hubbard model such 
a process takes place only via sequential tunneling, i.e. 
by a virtual disintegration of a pair. As a result, the cor- 
responding amplitude is suppressed with increasing U. 

In fact, a more probable mechanism for two-particle 
tunneling transitions takes place due to simultaneous 
tunneling of two particles. Yet, this type of process is not 
accounted for by the standard Hubbard Hamiltonian. In 
this Letter we show how such a direct two-electron tran- 
sitions can be included in the tunneling Hamiltonian by 
deriving the latter from the original Schrodinger equa- 
tion. We also demonstrate that such a modified Hubbard 
Hamiltonian produces an effective two-electron coupling 
due to on-site electron-electron repulsion. In this case 
pairs of electron with opposite spins can propagate along 
periodic systems, producing a current of entangled elec- 
trons. 

We start with an example of two electrons occupying 
a coupled-dot system, as shown schematically in Fig. 1. 
Each of the dots is represented by a square well potential, 
Vi,2{x), which contain a bound state at the same energy 
Eq. If both electrons occupy one of the dots (Fig. la,c) 
the total energy of this system is 2Eq + U. We assume 
that the inter-dot Coulomb repulsion is zero, so when- 
ever the electrons are in different dots (Fig. lb) the total 
energy becomes 2Eo. 

Such a system is usually described by the Hubbard 
tunneling Hamiltonian 

H = ^ {EoUis + ^Uisnig) - ^ r2o(aLa2s + H.c.) (1) 

i,s s 

where aj^ (a^s) creates (annihilates) an electron in the 
corresponding dot (i = 1,2) with the spin s (s — — s) 
and riis = al^ais- The amplitude Hq generates transitions 
between the dots via single electron tunneling. This am- 
plitude can be evaluated by using the Bardeen's formula 
as a product of two bound state wave functions at 
a point xq inside the barrier {a < xq < d — a): 



r!o = (k/to)$i(.to)$2(xo). (2) 

Here $1,2(2^) are the single electron wave functions in the 
left and the right dot: {K + Vi^2)^i,2 = £'o*^'i,2 with K is 
the kinetic energy operator, and k = (2m|i5o|)^/^. Note 
that the value of is weakly dependent on [Q . For 
instance, for square well potentials (Fig. 1) and a — > 0, 
the wave functions $i(a;) = k^^^ exp (— «;|a::|) and $2(2;) = 
K^/^ exp (— k|c? — x\), so that fig = {k^ /ra) exp(— kc?). 
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Fig. 1: Two electrons with opposite spins in a double- 
well potential V{x) = Vi{x) + V2{x). 

The wave function |^'(i)) describing a motion of two 
electrons inside the double-dot can be represented by 

E bu'{t)al^4i\0), (3) 

i,i' = l,2 

where &ii(<), 622 (t) are the probability amplitudes of find- 
ing two electrons in the first or in the second dot respec- 
tively (Fig. la,c). The probability amplitude 612 (i) corre- 
sponds to the configuration shown in Fig. lb, and &2i(0 
is the same but with an interchange of the spins. Substi- 
tuting Eq. (^) into the Schrodinger equation i9f |5'(t)) — 
H\'i>{t)), we obtain the following coupled equations for 
the amplitudes b{t) 
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ibii{22) = C^Eo + [/)6ii(22) + ^o{bi2 + 612) (4a) 



*^i2(2i) = 2i?o6i2(2i) + ^o{bii + 622) 



(4b) 



Let us solve these equations for the initial conditions 
&ii(0) = 1, 6i2(2i)(0) = 622(0) = 0, Fig. la. Con- 
sider first the case of no Coulomb repulsion between 
the electrons, U = 0. Then the probability of find- 
ing both electrons in the same dot, Pii{t) — \bii{t)\'^ 
(Fig. la,c), is Pii(i) = cos^{nat) and P22{t) = sin'*(f^oi)- 
Respectively, the probability of finding two electrons in 
different dots (Fig. lb) Pi2(t) = |&i2(0l^ + l^2i(i)P is 
PMt) =sm^i2not)/2. 

However, in the case of strong Coulomb repulsion, 
U ^ 4f2o, the time-behavior of two-electron system is 
very different. The solutions of Eqs. (jj) depend on a 
small parameter ^Iq/U. To second order of this parame- 
ter we find Pii{t) = cos'^{ujt), P22{t) = sin^(a;i) with uj ~ 
2VLI/U < rJo and Pi2(t) = {29.1/ U"^) sir? {U 1/2) < 1- 
Thus two electrons oscillate together between two wells, 
but with much lower frequency than in the noninteracting 
case. The probability of finding two electrons in different 
wells (Fig. lb) is strongly suppressed. 

The above result can be understood in the following 
way. If a dot is occupied by two electrons, the tunneling 
of one of these electrons out takes place with a variation 
of the total energy, /S.E = —U, Fig. 1. The probability of 
such a process is ~ fll/AE"^, and therefore it is strongly 
suppressed. This is analogous to the time-evolution of 
a single particle in a double-well potential. The particle 
remains in the same well where it was initially localized 
if the energy levels in the wells are misaligned. However 
if the levels are of the same energy, the particle oscil- 
lates between two wells (Rabi oscillations). In our case 
these are the two-electrons states. Fig. la,c, which have 
the same energy {2Eo + U). Therefore the two electrons 
oscillate between these states. Yet in the framework of 
the tunneling Hamiltonian Eq. (|l|) these oscillations are 
generated by a sequential tunneling (co-tunneling) via 
the intermediate state of Fig. lb, so that the frequency 
of oscillations is suppressed by a factor Hq/U. 

It follows from these arguments that the probability 
of two-electron tunneling between the states of Fig. la,c 
would be enhanced if both electrons tunnel simultane- 
ously, keeping the energy of the system constant. These 
processes must be accounted in any tunneling Hamilto- 
nian description by introducing the corresponding terms 
generating direct two-particle transitions. These terms 
are naturally arise whenever the Hubbard Hamiltonian is 
consistently derived from the Schrodinger equation. Such 
a derivation is presented below by using the Green's func- 
tion technique of the two-potential approach |^,^ . 

Consider the Schrodinger equation for two electrons in 
a double-well potential V{x) — Vi{x) + V2{x), Fig. 1. For 
simplicity we disregard the spin by treating two electrons 
with opposite spin as distinguishable particles. The total 
Hamiltonian can be written as H = Hi + H2 + Uc with 



Hi = K{xi) + V{xi) for i = 1,2, and Uc{xi — X2) is the 
Coulomb repulsive potential. We assume that Uc = U 
for |ri — r2 1 ^ c and drops down outside this region, where 
a <C c < d. 

We start with two electrons localized in the first dot, 
|*(0)) = where (xi, a;2|$ii) = $i(a;i)$i(a;2)- To 

determine the time-development of this system we apply 
the Laplace transform \^{E)) = \^{t)) eyi-p{iEt)dt. 
The Schrodinger equation then reads 



{E-H)\i!{E)) =i|*(0)) 



Using the Green's function technique 
the total wave function as 



(5) 



, we can rewrite 



\^{E)) = [1 + G{E)V2]bii{E)\^ii) (6) 
G{E)=Gi{E){l^kii)[l + V2G{E)], (7) 

where Gi{E) = [E-K-Vi)-^ with K = K{xi)+K{x2), 

y\,2 = Viaixi) + ^1^2(2^2) + C/c'^^(2;i - 2^2) and An = 
is the projection operator. Here we repre- 



sented Uc ~ + U(j , where Uq = Uc for xi,X2 ^ c 
and drops to zero outside this region. Correspondingly, 
Uc ior d — xi,d — X2 ^ c and drops to zero out- 



U. 



(2) 



side this region. The amplitude bii{E) in Eq. (|6|) is the 
Laplace transform of the probability amplitude feii(t) of 
finding both electrons in the first dot. Fig. la. One ob- 
tains from Eqs. (|)-(0) @ 



bii{E) = 



E-2E0-U- ($1111/2 + V2G{E)V2\<^ii) 



(8) 



Eq. (^) can be treated iteratively by expanding G, 
Eq. (^, in powers of Gi. Since E ~ 2Eo + U the Green's 
function Gi is dominated by the bound state pole in 
its spectral representation, Gi{E) |$ii)($ii|/(i^ — 
2Eo - U) for E ^ 2Eo + U. Yet, this pole is ex- 
cluded by the projection operator 1 — An. The remaining 
part of Gi gives rise to the corresponding (Born) series 
for \'^{E)) which, however, converges very slowly. We 
therefore look for a different expansion for G which con- 
verges much faster than the Born series. This can be 
achieved by expanding G in powers of the Green's func- 
tion G2{E) = 1/{E — K — V2) by using the relation 



G = G2(l + ViG) - G2Aii(l + V2G) . 



(9) 



Since the second dot contains the two-electron bound 
state at the same energy 2Eq -\- U &s the first dot, the 
Green's function G2{E) can be replaced by G2{E) 
1*22) ($22 1 /(£^ - 2£:o - U) for E ^ 2Eo + U, where 
(xi, 2:2 1 $22) — $2 (2^1) $2 (2:2 )• Substituting this expres- 
sion into Eq. (||) we find that in the contrast to Eq. (|^) 
the bound state pole is not cancelled by the projection 
operator 1 — An. In addition, the second term in Eq. (|^) 



2 



is proportional to the overlap (<&22|*I'ii) and can be con- 
sidered as a small perturbation. If we neglect this term 
Eq. (0) can be easily solved thus obtaining 



G{E) 



|$22)($2 



E-2E0-U - 2812 



(10) 



where ^12 — (<I'2|V^i|*i'2) is a small energy shift. Note 
that the Green's function G{E) contains also the poles at 
E ~ 2i?o, corresponding to separated electrons, Fig. lb. 
Since E ~ 2Eq + U these poles do not affect G{E) for 
large U. Substituting Eq. ^ into Eq. (|) we find 



\^{E)) = hii{E)\^ii) + b22{E)\^22) 



where 



h22{E) 



SI2 



E-2E0-U - 2<5i 



-buiE) 



(11) 



(12) 



is the Laplace transform of the probability amplitude 
b22it), Eq. (|h, and O2 = -(*22|V2|4'ii)- Then substi- 
tuting Eq. (pO|) into Eq. (H) and using Eq. (12) we obtain 



{E-2E0-U- 262i)bii{E) + n2b22{E) = I , (13) 

where 821 — ($i|V2|$i) (in our case 8x2 = ^21 ~ S). 

Performing the inverse Laplace transform, Eqs. (|l2|)- 
( |l3| ) can be rewritten in a matrix form as 



bn{t) 
622(0 



2E^ 



-Vt2 



U -^2 

2E', 



U 



bn{t) 

b22{t) 



(14) 



where Eq = Eq + S. Using the occupation number repre- 
sentation and re-introducing the spin variable one finds 
that Eq. (|l|) can be rewritten as the Schrodinger equa- 
tion 19(1* (t)} = iJ I *(*)), where |*(i)} is given by Eq. (|) 
and the Hamiltonian 

1=1,2 

- ft2{a\^a\^a2^a2i + H.c.) . (15) 

Here the amplitude il2 — — ($22|V2|*i'ii} generates a 
direct coupling in 2-dimensional Hilbert space between 
two-electron states. Fig. (la)^Fig. (Ic). This cou- 
pling can be calculated directly by using ($221^2 = 
($22|(2£'o — K) and integrating by parts. As a result we 
find the following simple expression (c.f. |||,0,^). 

^2 = — $i(a;o)$2(xo) / ^i{x')<i>2ix')dx' . (16) 
m J 

For the square-well potentials (Fig. 1) and a ^ 0, 
Eq. (|l|) yields fl2 = {fid/\Eo\)fll, where fla is single- 
electron tunneling amplitude, Eq. (|^). As expected 
the direct two-electron transitions, generated by fl2, 
dominate over the sequential two-electron transitions, 



^ 2^1q/U, with an increase of U. This dominance of 
course is not limited to square-well potentials. 

Note that Eq. (|l^) has been derived for long range re- 
pulsive potential Uc- For short range potential such that 
Uc{xi —X2) for \xi — X2I ^ a, the wave function $11 
and $22 are not given by a product of single electron 
wave functions. Yet the amplitude fl2 can be evaluated 
by using the Bardeen formula pjq ] . 

Comparing Eq. ( p^ ) with Eq. (|1|) wc find that the term 
^oo\s<i2s describing single electron transitions does not 
arise in Eq. ( p^ since the corresponding poles of G{E) 
at £■ ~ 2Eo were neglected in our derivation. However 
with a decrease of [/, such that SI2 ^ 2Qq/U, these poles 
must be taken into account. As a result we arrive to 
the Hamiltonian (^ supplemented with the last term of 
Eq. (15). Solving the Schrodinger equation with such 



a modified Hamiltonian for the initial conditions cor- 
responding to Fig. la we find in the limit of J7 ^ JIq 
that Pii(i) = cos2(u;'t), P22(t) = sin^(w'i) and Pi2(t) = 
{2nl/U^)sin^{Ut/2), where uj' = jf^a - {2nl/U)\. Note 
that the transition amplitudes flo and 1^2, given by 
Eqs. (H ) and (^6|) change the sign if the wave functions 
$1,2(2^) are of a different parity. 

Our procedure can be extended to a general case of 
N coupled wells. Detailed derivation will be presented 
elsewhere. Here we give our final result, representing the 
Hubbard-type Hamiltonian with additional terms for di- 
rect two-electron transitions 



N 



U 



N-l 



-^al^alsO-t+i-sai+i,s + H.c. 



(17) 



Here i = 1, . . . ,N and Eq = Eq + 6, where (5 is a small 
energy shift defined in the same way as in Eq. ([l^). The 
two-electron coupling Q2 is given by Eq. ( |l6|) for U :s> flo, 
but ^2-^0 for ?7 — > 0. In this limit, however, the contri- 
bution from direct two-electron transitions is suppressed 
with respect to that generated by single electron transi- 
tions. Thus Eq. ( [1^ ) with ^l2 given by Eq. ( [l^ ) can be 
used for any values of U. 

As in the previous case of a double-well potential, an 
electron pair with opposite spins initially localized in one 
of the wells cannot be separated in the limit of J7 — » 00 
Such a pair can only move as a whole object due to 
direct two-electron tunneling between neighboring wells 
(the amplitude of co-tunneling vanishes in this limit). As 
a result a mini-band of the width 2Q2 appears, providing 
a current of entangled electrons. 

If ^2 ~ 2^1^/11, the propagation of an electron pair is 
govern by the both co-tunneling and direct two-electron 
transitions. In order to assess their importance with re- 
spect to the single electron transport we solved numeri- 
cally the Schrodinger equation idt\'^{t))) = H\'i'{t)) with 
the Hamiltonian ( pT[ ) and the wave function given by 
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Eq. (|) with = I, N ioT N = 10 and U = 10% 
and the initial conditions corresponding to the first site 
occupied by two electrons. The results of our calculations 
are shown in Fig. 2. We plot there the probabilities of 
finding one and two electrons at the last site, P/v(<) and 
Pjvjv(i), for = O.ir^o and — 0. Although for the 
chosen values of parameters the value of ^2 = 0.1% is 
smaller than that of the co-tunneling one {2Uq/U), the 
direct two-electron transitions strongly affects the pair 
transport. For instance the pair reaches the last site con- 
siderably faster and with much larger probability. 



P ■ (t) 




Fig. 2: The probability of finding one and two electrons 
at the last site as a function of time. The solid lines cor- 
respond to % = O.lfio and the dotted lines to % = 0. 

It follows from Fig. 2 that the probability of finding a 
single electron at the last site is much lower than that 
of an entangled electron pair. Yet, one cannot estab- 
lish from this figure that this is true for all sites. More- 
over, it is natural to assume that the electron pair de- 
cays very fast. Indeed, the weight of all states with 
two electrons are at different sites is much larger then 
those occupying the same site. In order to investigate 
this point we evaluated numerically the total probabil- 
ity of finding two electrons at the same site inside the 
chain, Ppair{t) = X^i ■ The results for the pre- 

vious values of parameters are displayed in Fig. 3. It 
follows from this figure that contrary to the expectations 
the two-electron pair is stable during long time interval. 
It is interesting to note that the direct two-electron tran- 
sitions affect the stability of this pair very little. 

The entangled electron pair transport discussed above 
can be realized in multi-dot structures coupled to the 
reservoirs (emitter and collector) with a voltage bias 
larger than U . Then two electrons with opposite spins 
can enter simultaneously the first dot and then propagate 
across the entire structure. The single and the pair elec- 
tron currents in the collector can be found in the most 
simple way by using the Bloch-type rate equations for 
the reduced density matrix derived in |^ . 
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Fig. 3: Total probability of finding two electrons bound 
as a function of time, Ppair{t)- The solid and dotted lines 
correspond to % ~ 0.1% and % = 0, respectively. 

The two-electron coupling can also take place between 
electrons in neighboring dots via the inter-dot Coulomb 
repulsion. Indeed, if the separation energy of such a pair, 
Ui^i+i — Ui_i+2 is larger then a single-electron transition 
amplitude, %, the pair moves across the system as a 
whole object. This phenomenon for three coupled dots 
has been studied in but only taking into account the 
co-tunneling mechanism for the pair transport. 

In conclusion, we demonstrated that the Coulomb re- 
pulsion can bound electrons in multi-well periodic struc- 
tures, so that a current of entangled electrons with op- 
posite spins would appear. The propagation of such 
electron pairs across the system is dominated by direct 
two-electron transitions absent in the standard Hubbard 
model. We expect that these two-electron transitions 
play an important role in different physical processes. 

The author is indebted to M. Heiblum, Y. Imry, M. 
Kugler and A. Yacoby for fruitful discussions. 
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